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We study how the resonant decay of moduli fields arising in the Minimal Supersymmetric Standard 
Model (MSSM) could affect large scale curvature perturbations in the early universe. It has been 
known for some time that the presence of entropy perturbations in a multi-component system can 
act as seeds for the curvature perturbations on all scales. These entropy perturbations could be 
amplified exponentially if one of the moduli decays via stochastic resonance, affecting the curvature 
power spectrum in the process. By imposing the COBE normalization on this power spectrum, one 
could put constraints on the masses and couplings of the underlying particle physics model without 
having to rely on collider experiments. We discuss in detail the case of the MSSM but this method 
could be applied to other theories beyond the Standard Model. 



PACS numbers: 98.80.Cq 

I. INTRODUCTION 

Many models of physics beyond the Standard Model 
contain light scalar fields [H |2] which are usually called 
moduli fields. These field could have played an important 
role in the early Universe and therefore cosmology is a 
natural test-bed to study the properties of the moduli. 
Generically, models beyond the Standard Model predict 
the existence of multiple moduli fields. In such multi- 
component systems, one must consider the presence of 
entropy perturbations between the different constituents 
of the system in the early universe (see e.g. [3]). 

It is well-known that entropy modes (associated with 
non-adiabatic pressure perturbations) can act as a source 
for curvature perturbations on all scales ([4 and refer- 
ences therein) . Generally, the contribution of the entropy 
modes to the curvature power spectrum is small and can 
be neglected. However, under certain circumstances, the 
non-perturbative decay of one modulus field (e.g. the in- 
flaton itself) can exponentially amplify the amplitude of 
an entropy perturbation mode [H E] El E] • If this happens, 
the contribution to the curvature power spectrum arising 
from entropy perturbations could dominate over the adi- 
abatic modes (as in the curvaton scenario, see e.g. [B]). 
By demanding that the curvature power spectrum re- 
spect the bound set by the COBE normalization [9] , one 
can put constraints on moduli parameters (masses and 
couplings) , hence constraining certain models of physics 
beyond the Standard Model (see [10l E] for applications 
of this idea to string inflation models). 

The Minimal Supersymmetric Standard Model 



(MSSM) is a good example of a model having a large 
moduli space [12]. Indeed, the scalar potential of the 
MSSM has a large number of flat directions along which 
the potential vanishes |13j . Supergravity corrections usu- 
ally spoil the flatness of the potential [14 . Nevertheless, 
there exists a large class of models with non-minimal 
Kahlcr potential for which scalar fields can develop large 
vacuum expectation values (vevs) along flat directions 
during inflation |15| 1X6] - 

After inflation, the vev closely tracks the minimum 
of the potential as it evolves to lower scales |T7]. Once 
the Hubble parameter is of the order of the soft SUSY- 
breaking scale, the flat direction starts oscillating about 
its true minimum. Assuming an F-term interaction of the 
form AV = 9 \<p\ \x\ > where tp is the moving vev and 
X is some MSSM field, these oscillations could trigger 
the non-perturbative production of \ particles whenever 
the change in the frequency of \ ceases to be adiabatic 

KA4 > i) ng. 

However, since flat directions are inherently complex 
fields, they do not typically pass through |<^| = during 
their oscillations, ensuring that w x /w 2 <C 1 at all times 
[19] . Nevertheless, under certain circumstances, the ellip- 
ticity of the trajectory in field space can be small enough 
such that one can apply the well-know results from para- 
metric resonance [20] . 

Furthermore, it was recently shown [2TJ[22 I 23, 24 that 
if two or more flat directions are simultaneously excited, 
the rotation of the vevs could trigger exponential produc- 
tion of MSSM fields. Thus, even if the ellipticity of the 
trajectory is large, we expect entropy modes between the 
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flat direction condensates and the created particles to be 
amplified. 

In this paper, we focus on the simpler case where the 
motion of the oscillating vcv is purely radial. We con- 
sider the case of a quantum field x whose mass during 
inflation is somewhat smaller than the Hubble scale and 
which is coupled to a flat direction. During inflation, 
quantum fluctuations of this field are stretched to super- 
Hubble scales and form an almost scale-invariant spec- 
trum at the end of inflation. The subsequent decay of 
the flat direction triggers an exponential amplification of 
these super-Hubble modes via stochastic resonance [15] , 
These modes correspond to an entropy perturbation that 
then acts as a source for curvature perturbations. We 
compute the curvature power spectrum and compare it 
to the bound set by the COBE normalization. Note that 
induced curvature fluctuations from entropy modes in su- 
persymmetric models were recently also studied in |25j . 
New to our work is the focus on the parametric amplifica- 
tion of the entropy mode during the phase of oscillation 
of a flat direction field. 



II. THE MODEL 

We consider a MSSM flat direction parametrized by 
the order parameter tp. Soft SUSY-breaking and Hubble- 
induced mass terms as well as superpotential corrections 
of the type W D Xtp n /nM n ~ 3 (n > 3) lift the flat direc- 
tion by giving it a potential of the form [17 : 



V(tp) 



(ml-c H H 2 M 2 - 



2 M 2n ~ 2 

l A ! 2 M 2«-6' W 



where M is the scale at which the non-renormalizable 
terms become important. Here we neglected soft SUSY- 
breaking and Hubble-induced A-terms. With a minimal 
choice of Kahler potential, ch — —3 and the flat direction 
mass is positive. Therefore, tp settles at the origin during 
inflation and does not lead to any interesting cosmologi- 
cal consequences. Nevertheless, there exists a large class 
of Kahler potentials for which cu > [THl HE1 [SB] . In 
these cases, the potential (JT|) admits a minimum at: 



2 s J{n-l)c H HM;~ 
|A| 



(2) 



where we took the cutoff to be the Planck scale, M = M p . 
During inflation, the flat direction settles down at this 
minimum. The renormalizable MSSM superpotential |13j 
naturally provides couplings between the modulus field 
tp and other MSSM fields Xi- The generic form of the 
potential for a field coupled to flat direction via an F- 
term interaction is: 



V(x) = \™lx 2 + \g 2 W\ 2 x 2 . 



(3) 



In the following, we will always assume that m x <C g\ip\ 
such that we can neglect the mass of the x field. 



During inflation, we demand that g 2 \tp m i n \ 2 < Hj in 
order for x to acquire un-suppressed fluctuations on large 
scales. This can be achieved by choosing the smallest 
possible value of n, namely n = 4. Indeed, |y| m j n lies 
parametrically between H and M p and it is easy to see 
that for large n, \tp\ m i n —> M p . A way to relax this con- 
dition would be to consider \J2cnj\ <C 1. However, we 
need cu ~ [17] in order for the quantum fluctua- 

tions of x n °t to dominate the mass of tp (otherwise tp 
would not oscillate freely). This in turn implies A ^> 1 
in order to have \J2cnj\ -C 1 which can be regarded as 
a fine-tuning. Hence, for the rest of this paper, we focus 
on flat directions that are lifted at the n — 4 level (i.e. 
the smallest possible value of n). For this particular case, 



the condition g 2 \tp m in\ 2 < Hf translates to 
g 2 < Hj/Mp . 



(4) 



where Hi is the Hubble parameter during inflation. 

After inflation, the flat direction closely tracks the po- 
sition of the minimum as it shifts to lower scales be- 
cause of the time-dependence of the Hubble parameter. 
Once chH 2 ~ m 2 , the flat direction gets a positive mass 
squared and starts oscillating with frequency m v . Since 
we are taking A-terms to be negligible, these oscillations 
are for all practical purposes one dimensional. They lead 
to a time-dependence of the x field mass that can trigger 
non-perturbative amplification of x fluctuations when- 
ever the change of the x frequency u> x — g\tp\ ceases to 
be adiabatic [18] . We take the universe to be dominated 
by radiation after inflation. In a radiation dominated 
universe, the time evolution of the flat direction once it 
starts oscillating is given by: 



¥>(*) 



(/?o sin (m v t) _ 
(rV) 3 / 4 



= $(t) sin (m v t) t 



(5) 



where tp is given by ^ evaluated at H(t) ~ m v . 

On the other hand, the perturbations in the x field 
obey the wave equation in a FPW universe with metric 
ds 2 =dt 2 -a 2 (t)dx 2 : 



Xk + ZHxk+[%+g 2 W)\ z )xk = 



(6) 



where a is the usual scale factor. We can simplify this 
equation by rescaling the field as X^{t) = a 3 / 2 (t)xk(t) 
|18j . This leads to the much simpler equation: 



where 



a 2 (t) 



X k + uj{X k = 0, 



+ g 2 \^t)\ 2 --H 2 - 



3 a 
2a' 



(7) 



(8) 



In a radiation dominated universe, the two last terms of 
(fsp add up to +3iJ 2 /4 and thus ui\ is strictly positive. 
The oscillations of the flat direction may trigger the am- 
plification of quantum fluctuations in x y i a stochastic 
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resonance. Such a resonance will occur if the condition 

^H9 2 M p 



2 2 

9 Po 
4m2 



2|A|7 



» 1 



(9) 



is satisfied. Since we expect c H ~ C(l) (see above) and 
taking |A| ~ 0(1) (which yields a lower bound on q), this 
condition is realized provided that 



g 2 M p > m ¥ 



(10) 



at the beginning of the resonance. Assuming that this 
condition holds, the solution to ^ is given by [TB] : 



\Xk(t)\ 



\Xo(k)\e" km ^ 



(11) 



where Xo(&) is the initial amplitude of the fc*' 1 mode of 
the x fluctuations at the beginning of the resonance and 
fik is the effective Floquet exponent (whose value is dis- 
cussed later). This growth corresponds to the squeezing 
of the quantum vacuum state of the x field. The squeezed 
vacuum state is highly excited and can be viewed as con- 
taining a condensate of x particles with momentum k. 
The solution ( 11 ) is only valid for modes satisfying: 



< 



gm v <S>(t) 



(12) 



For modes outside this resonance band, the adiabaticity 
condition is not violated. Hence, the effective Floquet 
exponent vanishes and no growth occurs. Note that large 
scale modes are preferably amplified compared to short 
wavelength modes. 

The resonant amplification takes place until backreac- 
tion shuts off the resonance. This occurs when the mass 
of the oscillating flat direction begins to be dominated by 
the fluctuations of the x held, that is, when g 2 (x 2 ) — m %- 

To summarize the results of this section, we have shown 
that the oscillating homogeneous tp field leads to stochas- 
tic resonance of the x ^sld in the same way that the 
background homogeneous inflation field leads to stochas- 
tic resonance of matter fields which it couples to [T8|I20|. 
generating a squeezed state of low frequency x particles. 
The resonance of x eventually shuts off the resonance by 
changing the time evolution of the background tp field. 



III. CHANGE IN THE CURVATURE 
PERTURBATION ON LARGE SCALES 

We now turn our attention to the effect of the am- 
plification of the x fluctuations on large scale curvature 
perturbations. In the presence of a non-adiabatic pres- 
sure perturbation Jp na d, the curvature perturbation on 
uniform-density hypersurfaces £ varies according to the 
equation [271IS5] : 



C 



H 

P + P 



8p, 



nad ; 



(13) 



where the non-adiabatic pressure perturbation is defined 
as: 



Sp. 



nad 



Sp — -Sp. 

P 



(14) 



In our case, the entropy mode is generated by the field 
X and is quadratic in x (there is no rolling homogeneous 
X background). The leading-order pressure and energy 
density perturbations are given by: 



c 1 2 2 2 , 1 -2 , c 

Sp = 2# V * + 2 X + Pa 

s 2 2 2 , 1 -2 , r 

°P = -~9 <P X + 2* +°Pa, 



(15) 
(16) 



where we defined tp = \<p\ (we are considering radial mo- 
tion of the field only) . Sp a and Sp a are the adiabatic com- 
ponent of the perturbations, that is, Sp a — {pol Po)5p a - 
The adiabatic component of the fluid is made up of the 
field tp and of the radiation produced during inflationary 
reheating. Thus, the background density and pressure 
are given by: 



p = (p 2 + m 2 <p 2 + p Q 



• 2 2 2 

p = tp - m v (p 



Po, 



(17) 
(18) 



where po and po are the background pressure and energy 
density, subtracting the contribution from tp. Since we 
are assuming that the flat direction decays during the 
post-inflationary phase of radiation domination, we have 
Po = (l/3)jOo- Combining the above results, we obtain 
the following expression for the non-adiabatic pressure 
perturbation: 



Spnad = 



tp(g 2 X 2 ^ - + (^) {\X 2 - ffW) 

3Htp 2 + 2Hp 



2tp5p a 



1m%tp 



'■P 



3 3Htp 2 + 2Hp Q 



(19) 



Since we are interested in the growth of the non-adiabatic 
pressure perturbation at the time when the modulus 
field starts oscillating, we have H ~ m v such that 



pa ~ m^M 2 . Meanwhile, using |2j), the time-averaged 
time derivative of the modulus field is {tp 2 ) ~ m? M p . 
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Therefore, 



(v 2 



Pa 



M„ 



« 1, 



(20) 



and we can safely neglect the first term in the denomina- 
tor of (19). Similarly, the second term of the numerator 



dominates over the first one since it is proportional to 
po- Moreover, the last term of (19 1 is independent of \ 



(the field whose variance is growing exponentially) and 
we thus expect it to give a negligible contribution to the 
non-adiabatic pressure. Keeping only the dominant part 
to the non-adiabatic pressure perturbation, the equation 
of motion for the curvature perturbation reads: 



c 



1 



1 



6HM 2y 2 



I ■ A 2 2 2\ 



where we used p = 3M 2 H 2 . Using \x\ 2 ^ g 2 <P 2 \x\ 2 an d 
H = l/2t, this equation can readily be integrated: 



Cnad 



g 2 2 f Xf sin 2 a; 

'Xo| 



6M p m v 



e 2fl! * x dx, (22) 



where x — m v t and where we used ^ and (111 to 
obtain the time dependence of the fields. The integral 
runs over the duration of the resonance and the value of 
x f corresponds to the time when the backreaction of the 
produced \ particles on the oscillating modulus becomes 
important, that is, when g 2 (x 2 ( x f)) — TO ^>- Given the 
value of the Floquet exponent and the duration of the 
resonance, the integral in ( |22[ ) can be done numeri- 
cally and yields a pure number that we shall denote by /. 

The power spectrum of the non-adiabatic curvature 
perturbations is given by: 



V, 



Cnad 



:\(nad(k)\ 2 , 



(23) 



where Q n ad{k) is related to Q na d by a Fourier transform. 
Using ( 22 1 and performing the Fourier transform integral 
yields: 



k 3 ( 



2 I 



P^(\k'\)V^(\k-k'\) 
Ik'Plk-k'F 



-}:: \ !)_U „/// - / ./ IV | • .IV — .IV 

(24) 

where "P£ es is the power spectrum of the \ particles at 
the beginning of the resonance which is given by: 



1 X 



El 

2tt 



Hi 



k 



2,1 



where k en( i is the last comoving scale to have exited the 
Hubble radius at the end of inflation and 5 = g 2 f 2 /3iff , 
where the / subscript refers to a quantity evaluated dur- 
ing the inflationary era. The first term is term is the 
usual scale invariant spectrum for a field that is light 
during inflation, the second one accounts for the Hubble 



damping between the end of inflation and the beginning 
of the resonance while the k-dependent term takes into 
account the small mass of the \ field during inflation. 
Since we have 6 <C 1, we can neglect the slight scale 
dependence of the power spectrum and perform the mo- 
mentum integral. This integral has a IR divergence that 
must be regulated. A natural infrared cutoff is provided 
by the first comoving mode that exited the Hubble scale 
at the beginning of inflation. Taking this mode to cor- 
respond to our present Hubble radius, we are basically 
integrating over today's visible universe, i.e. the only 
patch of the universe that we have information about. 
The resonance band (12 1 naturally provides a UV cut- 



off for the non-adiabatic power spectrum (although the 
answer depends very mildly on this choice since the in- 



(21) tegrand in (24 1 is rapidly converging). The answer then 



depends logarithmically on the total number of e-foldings 
of inflation. Taking N = 60 and imposing the constraint 
g 2 < Hi/Mp, we obtain an upper bound on the non- 
adiabatic power spectrum: 



TO 2 / g 2 I 

2ir 2 \&M V 



Hi 
2tt 



< 



Hi 



(26) 



Given an inflationary scale Hi , the overall size of the non- 
adiabatic power spectrum is thus determined by the en- 
hancement factor due to the resonant decay of the modu- 
lus field. Therefore, in order to see if this power spectrum 
could exceed the bound set by the COBE normalization, 
we need to obtain as estimate of the value of / which in 
turn depends on the length of the resonance and on the 
value of the effective Floquet exponent. 



IV. STRENGTH OF THE RESONANCE 

To estimate the resonance enhancement factor, one 
needs to determine the value of the effective Floquet ex- 
ponent. On large scales, we can neglect the scale depen- 
dence of fj,k, and its approximate value is then given by 
UHl: 



1 

27 



In 3 ~ 0.17. 



(27) 



The length of the resonance is obtained by solving the 
equation g 2 (x 2 ( x f)) — m % f° r x f- The variance for the 



X field at the end of the resonance is 

3 Hf m w ( a(x ) x 3 " 



(25) (x 2 (x f )) 



,2n(xf-x ) 



(28) 



Stt^V? Hj \a(x f ). 

Using Q to calculate the amplitude if during inflation 
and noting that a(t) oc i 1 / 2 , we can obtain an approxi- 
mate solution (up to a logarithmic correction) for xf. 



-—In 
2[i 



El 

2ir 



+ 0(lnx f ) (29) 
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which is valid if xj ^> 1, that is, if m v ^> Hj/M p , which 
is also the condition that must be satisfied to have a long- 
lasting resonance. Note that this result is independent of 
the coupling constant g. 

Combining this result with the requirement that the 
mass of tp is negligible during the period of inflation and 
with the constraints on g found earlier (see |4]and 10 ), we 
obtain the following parameter space for which resonance 
could possibly lead to a curvature power spectrum that 
exceeds current observational bounds: 



H 2 

Hi > m v > —f- 



Ml 



< 9 < 



(30) 
(31) 



The crucial question now is: given values of parameters 
lying within the above bounds, is it possible for the res- 
onance to lead to an observable effect on the curvature 
power spectrum? First, to obtain an upper bound on 
the duration of the resonance, we take m v to be close to 
the inflation scale, that is, we write m v = eHj. We can 



now integrate ( 22 ) and calculate the upper bound on the 



non-adiabatic curvature power spectrum as a function of 
the inflationary scale. The results are shown in Figure [T] 
for various values of e. 

We see that the non-adiabatic power spectrum can 
exceed the limit set by COBE [9] only if the mass of 
the modulus field is close to the inflationary scale and if 
the latter lies close to the Planck scale. For example, if 
m v ~ 10 -3 Hj, then the non-adiabatic power spectrum 
could exceed observational bound only if Hj > 10 _2 M p . 
Lower moduli masses or inflationary scales yield non- 
adiabatic power spectra that do not contradict observa- 
tional bound and therefore do not allow us to put con- 
straints on the particle physics parameters. For moduli 



masses larger than 10~ 3 Hi, there is a region of the (m, g) 
parameter space which is excluded, provided that the in- 
flationary scale is larger than 10 _4 M p . For example, for 
Hi ~ 10 _3 M p the following values could be excluded: 



io- 3 > » io- 6 



Mr, 



M r , 



1/2 



< g < !Q- 3/2 . 



(32) 
(33) 



In particular, the usual value ra v ~ O(TeV), which 
allows SUSY to solve the hierarchy problem, is permitted 
for all values of the inflationary scale. 



V. DISCUSSION AND CONCLUSION 

We have shown that for a MSSM modulus field which 
acquires a negative Hubble-induced mass during inflation 
and subsequently decays via broad resonance, the growth 
of the curvature perturbations due to the non-adiabatic 
pressure does not contradict current observational bound 
for most of the parameter space. The non-perturbative 
decay of the modulus field could yield an observable effect 
only if the scale of inflation is very close to the Planck 
scale and if the soft SUSY-breaking mass is close to the 
inflationary scale. Consequently, the technique devel- 
oped above does not allow one to put constraints on the 
MSSM at the TeV scale. However, this result depends 
closely on the shape of the modulus potential . In par- 
ticular, the effect can only be large if the energy stored in 
the moduli field is a considerable fraction of the overall 
energy density of the universe. 



In this paper, we treated (inaccurately) the modulus 
field as a real field. However, MSSM flat directions are 
intrinsically complex fields and one could ask if our result 
apply if we consider the complex nature of the moduli. 
In references [2TJ [52 [331 IMj , it was shown that by taking 
into account their complex nature, flat directions can de- 
cay non-perturbatively if at least two flat directions are 
excited at the same time. Therefore, it seems reasonable 
to expect that parametric amplification of entropy modes 
(and by consequence of the curvature perturbations) can 
happen for more realistic model parameter values if we 
were to consider the full complex nature of flat directions. 

In conclusion, we point out that the large scale cur- 
vature perturbations that were amplified by stochastic 
resonance could be the dominant contribution to the cur- 
vature power spectrum. Indeed, within the model devel- 
oped above, inflation does not have to generate pertur- 



bations of the right order of magnitude. Instead, these 
could be generated at later times via the stochastic decay 
of a modulus field. 
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Upper Bound on Non-Adiabatic Power Spectrum versus Inflationary Scale 
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FIG. 1: Upper bound for the non-adiabatic power spectrum 
for different values of the modulus mass. For m v ~ 10 _1 Hi, 
the resonance could lead to an observable effect only if Hi > 
10~ 4 M„ while for m v ~ 10~ A Hi, one needs Hi > 10" 2 M„. 
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